本发明公开了一种基于海底数字水深模型特征提取的多波束水深图构建方法。方法包括基于原始多波束测深的DDM的构建;基于DDM建立坡度与二阶导数复合模型;基于DDM分块以提取特征点;基于特征点查询原始多波束水深;基于图层综合制作多波束水深图。本发明采用分块方式,基于DDM快速提取能表征海底地形基本特征的水深点,进而基于特征点坐标设置查询半径、在原始水深数据中快速查询与之最接近的点重新建模,克服了原始多波束水深数据量大难以提取特征的缺陷,也解决了直接使用DDM水深点或中央波束点难以准确反应原始探测水深的缺点。该方法在海底地形地貌制图、海洋测绘、海洋信息系统、计算机图形学与海底科学研究中具有非常重要的实际应用价值。
1.一种基于海底数字水深模型特征提取的多波束水深图构建方法,其特征在于,包括下列步骤:步骤1、基于原始多波束测深的数字水深模型DDM的构建,使用多波束后处理系统,对多波束原始数据进行如下处理:
1.1)对获取的原始多波束测深数据集合Raw={raw_i}_i=1,n进行潮位改正、吃水改正和声速改正,以消除水位和声速对多波束测深MBES的影响,形成初步处理后的数据集合Proc1={proc1_i}_i=1,n,n为自然数;
1.2)建立地形趋势基面,采用自动滤波方法对Proc1={proc1_i}_i=1,n进行处理,消除粗差对MBES的影响,形成后处理数据集合Proc2={proc2_i}_i=1,n;
1.3)采用三视图投影方式,处理Proc2={proc2_i}_i=1,n,进一步消除多种来源海洋噪声对MBES的影响,进而形成全波束离散的水深数据Proc3={x_i,y_i,z_i}_i=1,m,m为自然数;
1.4)采用距离反比加权方法IDW(Inversexa0Distancexa0Weighted)方法,对Proc3={x_i,y_i,z_i}_i=1,m进行数字建模,形成数字水深模型 DDM = { ddm k , l = { x k , l , y k , l , dep k , l } } k = 1 , K l = 1 , L , ]]>K为该模型行数,L为该模型列数,K和L均为自然数,x_k,l,y_k,l和dep_k,l分别为模型数据点的坐标和水深值;步骤2、基于DDM建立坡度与二阶导数复合模型,
2.1)采用八邻点方法,计算每个模型网格点的坡度值slp_(i,j): slp ( i , j ) = 1 8 Σ i = k k + 1 Σ j = 1 l + 1 a tan ( Δz Δd ) ]]>Δz=|dep_(i,j)-dep_(k,l)|,Δz为模型中编号为(i,j)和(k,l)两点的水深差值;ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="700"/>Δd为模型中编号为(i,j)和(k,l)两点间的距离;
2.2)采用八邻点方法,计算每个模型点的二阶导数值sec_(i,j): sec ( i , j ) = 1 8 Σ i = k k + 1 Σ j = 1 l + 1 a tan ( Δδ Δd ) ]]>Δδ=|slp_(i,j)-slp_(k,l)|
2.3)形成复合数字模型 DSSM = { dss i , j = { ddm i , j , slp i , j , sec i , j } } i = 1 , K j = 1 , L ; ]]>步骤3、基于DDM分块以提取特征点,
3.1)模型分块:以边长为d的正方形隔行错位对模型DSSM分块;d根据制图比例尺确定d=[0.018×scale,0.03×scale],其中scale为比例尺值;或者d由用户直接指定;隔行错位是指第i行以d为间距分块模型网格,至第i+Δi行时,首先空隔0.5d倍的间距,然后以d为间距分块模型网格;隔行错位,形成菱形分块效果;通过分块,获得新的菱形网格模型 DSSM 1 = { dss 1 I , J = { dss i , j } i = I , I + Δi j = J , J + Δi } I + 1 , K 1 J = 1 , L 1 , ]]>每个分块dss1_I,J由子模型网格ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="275"/>组成,K1和L1为分块模型的行数和列数,均为自然数;
3.2)分块模型数值计算:
3.2.1)计算 dss 1 I , J = { dss i,j } i = I , I + Δi j = J , J + Δi ]]>的平均水深 dep ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi dep ( i , j ) ]]>、平均坡度值 slp ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi slp ( i , j ) ]]>和平均二阶导数值 sec ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi sec ( i , j ) , ]]>Num为每个分块模型中的格网点数,是自然数;
3.2.2)按照距离反比加权方法IDW获取每个分块模型的中心点水深值Δdep_I,J;
3.3)分块模型凹凸性判断,
3.3.1)当ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="332"/>且ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="198"/>时,该分块模型曲面是下凹的;
3.3.2)当ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="330"/>且ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="208"/>时,该分块模型曲面是上凸的;
3.3.3)其他情况下,该分块模型难以确定凹凸性;
3.4)分块模型特征水深确定,
3.4.1)当分块模型是凹曲面时,遍历模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="494"/>获取最大水深点dss_{i_max,j_max};如有多个最大水深点时,选取坡度值slp_i,j最大且靠近中心的水深点;
3.4.2)当分块模型是凸曲面时,遍历模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="494"/>获取最小水深点dss_{i_min,j_min};如有多个最小水深点时,选取坡度值slp_i,j最大且靠近中心的水深点;
3.4.3)其他情况时,选取模型中心点dss_{i_cen,j_cen}作为特征点;
3.5)模型特征点标识,当dss_i,j为特征水深点时,设置其标记为1,否则设置其标记为0;步骤4基于特征点查询原始多波束水深,
4.1)遍历数据集合Proc3={x_i,y_i,z_i}_i=1,m,按照每个水深点的坐标值(x_i,y_i)查询它作用的模型网格点dss_I,J;模型网格点dss_I,J的坐标为(X_I,Y_J);
4.2)当 dis = ( x i - X I ) 2 + ( y j - Y J ) 2 < 0.5 × grid _ d ]]>时,且当dss_I,J为特征标识值为1时,使用IDW方法,基于数据集合Proc3={x_i,y_i,z_i}_i=1,m,重新计算模型点的水深值dep_new_I,J;grid_d为模型最近两点的空间间距;步骤5基于图层综合制作多波束水深图,使用成图系统或地理信息系统,按照等深线层、水深层和修饰层三个分层来制作来多波束水深图;
5.1)在等深线层中,根据成图比例尺,按照国标要求构建首曲线、辅助线、记曲线和记曲线上的水深标注值;
5.2)在修饰层中,添加图名、比例尺、图例、图框、经纬网线和海岸线;
5.3)在水深层中,遍历模型DSSM,当模型点dss_I,J特征标识为1时,且当水深值dep_new_I,J不为空值时,使用该点为特征点,当dep_new_I,J为空值时,使用dep_I,J作为特征点;
5.4)当最水深图中水深点过于稀疏时,减小模型分块大小值d,返回步骤3;
5.5)当最水深图中水深点过于密集时,增加模型分块大小值d,返回步骤3;
5.6)多图层叠加后,如出现图形要素相互叠加而不可见现象时,或者图形要素太密集时时,则删除一些辅助的图形要素。
2.根据权利要求1所述的基于海底数字水深模型特征提取的多波束水深图构建方法,其特征在于,步骤5.4)所述的稀疏为所制作的水深图中任意一点与其相邻点的距离D>2d;步骤5.5)和5.6)所述的密集为制作的水深图中任意一点与其相邻点的距离D<0.5d。技术领域
本发明涉及海底地形地貌制图、海洋测绘、海洋地理信息系统、计算机图形学和海底科学等技术领域。
背景技术
水深注记在海底地形制图、海图制作或海洋工程制图上是一项重要内容,通过水深图能定量展示海底地形、地貌特征,能直观地反应海区的通航能力,无论是在科学研究还是在实际航海中水深图都起着其他图件无法替代的作用。水深图与一般海底地形图的区别在于:以水深注记为主,等深线为辅,结合底质注记加以表示,水深是反映海底地势起伏的主要方法。传统的单波束回声测深获取的水深数据较少,在制作水深图时尽量保留大部分原始测深数据,其显著缺点是数据稀疏,导致海底地形特征点常被遗漏。目前在高要求的海底探测中已很少使用单波束测深。多波束测深MBES(Multi-Beamxa0Echoxa0Sounding)是当前最先进的水深测量技术,通过该技术获取的是海底全覆盖、高精度水深数据,基于该种数据可以制作大比例尺水深图。但多波束测深的显著特点就是海底探测的水深分辨率高、数据量大,因此在制作水深图时必须从海量测深数据中按照水深注记原则挑选一批符合制图要求的原始水深特征数据。目前基于多波束测深数据制作水深图的方法主要有两种:直接使用格网数据编制;使用中央波束编制。这两种方案均存在局限。直接用格网数据编制的缺点是没有采用实测数据,导致图件的可信度降低;而基于中央波束编制水深图的易遗漏测线间需注记的水深特征点。如何充分使用原始的多波束探测水深数据来制作多波束水深图目前尚未查询到科学的解决方案或文献。
发明内容
本发明针对存在多波束测深数据量大难以准确制造水深图问题,提出一种基于海底数字水深模型DDM(Digitalxa0Depthxa0Model)特征提取的多波束水深图构建方法,采取如下步骤构建多波束水深图:1、基于原始多波束测深的DDM的构建;2、基于DDM建立坡度与二阶导数复合模型;3、基于DDM分块以提取特征点;4、基于特征点查询原始多波束水深;5、基于图层综合制作多波束水深图。本发明是通过下述技术方案得以实现的:一种基于海底数字水深模型特征提取的多波束水深图构建方法,包括下列步骤:步骤1、基于原始多波束测深的数字水深模型DDM的构建使用多波束后处理系统,对多波束原始数据进行如下处理:1.1)对获取的原始多波束测深数据集合Raw={raw_i}_i=1,n进行潮位改正、吃水改正和声速改正,以消除水位和声速对多波束测深MBES的影响,形成初步处理后的数据集合Proc1={proc1_i}_i=1,n,n为自然数;1.2)建立地形趋势基面,采用自动滤波方法对Proc1={proc1_i}_i=1,n进行处理,消除粗差对MBES的影响,形成后处理数据集合Proc2={proc2_i}_i=1,n;1.3)采用三视图投影方式,处理Proc2={proc2_i}_i=1,n,进一步消除多种来源海洋噪声对MBES的影响,进而形成全波束离散的水深数据Proc3={x_i,y_i,z_i}_i=1,m,m为自然数;1.4)采用距离反比加权方法IDW(Inversexa0Distancexa0Weighted)方法,对Proc3={x_i,y_i,z_i}_i=1,m进行数字建模,形成数字水深模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="700"/>K为该模型行数,L为该模型列数,K和L均为自然数,x_k,l,y_k,l和dep_k,l分别为模型数据点的坐标和水深值;步骤2、基于DDM建立坡度与二阶导数复合模型2.1)采用八邻点方法,计算每个模型网格点的坡度值slp_(i,j): slp ( i , j ) = 1 8 Σ i = k k + 1 Σ j = l l + 1 a tan ( Δz Δd ) ]]>Δz=|dep_(i,j)-dep_(k,l)|,Δz为模型中编号为(i,j)和(k,l)两点的水深差值;ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="700"/>Δd为模型中编号为(i,j)和(k,l)两点间的距离;2.2)采用八邻点方法,计算每个模型点的二阶导数值sec_(i,j): sec ( i , j ) = 1 8 Σ i = k k + 1 Σ j = l l + 1 a tan ( Δδ Δd ) ]]>Δδ=|slp_(i,j)-slp_(k,l)|2.3)形成复合数字模型 DSSM = { dss i , j = ( ddm i , j , slp i , j , sec i , j ) } i = 1 , k j = 1 , L ; ]]>步骤3、基于DDM分块以提取特征点3.1)模型分块:以边长为d的正方形隔行错位对模型DSSM分块;d根据制图比例尺确定d=[0.018×scale,0.03×scale],其中scale为比例尺值;或者d由用户直接指定;隔行错位是指第i行以d为间距分块模型网格,至第i+Δi行时,首先空隔0.5d倍的间距,然后以d为间距分块模型网格;隔行错位,形成菱形分块效果;通过分块,获得新的菱形网格模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="700"/>每个分块dss1_I,J由子模型网格ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="270"/>组成,K1和L1为分块模型的行数和列数,均为自然数;3.2)分块模型数值计算:3.2.1)计算 dss 1 I , J = { dss i , j } i = I , I + Δi j = J , J + Δi ]]>的平均水深 dep ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi dep ( i , j ) , ]]>平均坡度值 slp ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi slp ( i , j ) ]]>和平均二阶导数值 sec ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi sec ( i , j ) , ]]>Num为每个分块模型中的格网点数,是自然数;3.2.2)按照距离反比加权方法IDW获取每个分块模型的中心点水深值Δdep_I,J;3.3)分块模型凹凸性判断3.3.1)当ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="330"/>且ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="201"/>时,该分块模型曲面是下凹的;3.3.2)当ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="336"/>且ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="200"/>时,该分块模型曲面是上凸的;3.3.3)其他情况下,该分块模型难以简单确定凹凸性;3.4)分块模型特征水深确定3.4.1)当分块模型是凹曲面时,遍历模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="498"/>获取最大水深点dss_{i_max,j_max};如有多个最大水深点时,选取坡度值slp_i,j最大且靠近中心的水深点;3.4.2)当分块模型是凸曲面时,遍历模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="488"/>获取最小水深点dss_{i_min,j_min};如有多个最小水深点时,选取坡度值slp_i,j最大且靠近中心的水深点;3.4.3)其他情况时,选取模型中心点dss_{i_cen,j_cen}作为特征点;3.5)模型特征点标识当dss_i,j为特征水深点时,设置其标记为1,否则设置其标记为0;步骤4基于特征点查询原始多波束水深4.1)遍历数据集合Proc3={x_i,y_i,z_i}_i=1,m,按照每个水深点的坐标值(x_i,y_i)查询它作用的模型网格点dss_I,J;模型网格点dss_I,J的坐标为(X_I,Y_J);4.2)当 dis = ( x i - X I ) 2 + ( y j - Y J ) 2 < 0.5 × grid _ d ]]>时,且当dss_I,J为特征标识值为1时,使用IDW方法,基于数据集合Proc3={x_i,y_i,z_i}_i=1,m,重新计算模型点的水深值dep_new_I,J;grid_d为模型最近两点的空间间距,使用0.5×grid_d重新建模目的在于使重建模型点水深接近真实值,不直接使用原位点探测水深是为了避免偶然探测误差对制图的影响;步骤5基于图层综合制作多波束水深图使用成图系统或地理信息系统,按照等深线层、水深层和修饰层三个分层来制作来多波束水深图;5.1)在等深线层中,根据成图比例尺,按照国标要求构建首曲线、辅助线、记曲线和记曲线上的水深标注值;5.2)在修饰层中,添加图名、比例尺、图例、图框、经纬网线和海岸线;5.3)在水深层中,遍历模型DSSM,当模型点dss_I,J特征标识为1时,且当水深值dep_new_I,J不为空值时,使用该点为特征点,当dep_new_I,J为空值时,使用dep_I,J作为特征点;5.4)当最水深图中水深点过于稀疏时,减小模型分块大小值d,返回步骤3;5.5)当最水深图中水深点过于密集时,增加模型分块大小值d,返回步骤3;5.6)多图层叠加后,如出现图形要素相互叠加而不可见现象时,或者图形要素太密集时时,可删除一些辅助的图形要素。所述的基于海底数字水深模型特征提取的多波束水深图构建方法,步骤5.4)所述的稀疏为所制作的水深图中任意一点与其相邻点的距离D>2d;步骤5.5)和步骤5.6)所述的密集为制作的水深图中任意一点与其相邻点的距离D<0.5d。本发明的有益效果本发明公开了一种基于海底数字水深模型特征提取的多波束水深图构建方法。针对目前多波束探测中数据量大、难以使用原始水深制作水深图的问题。该技术方法采用分块方式,基于DDM快速提取能表征海底地形基本特征的水深点,进而基于特征点坐标设置查询半径、在原始水深数据中快速查询与之最接近的点重新建模,克服了原始多波束水深数据量大难以提取特征的缺陷,也解决了直接使用DDM水深点或中央波束点难以准确反应原始探测水深的缺点。该方法在海底地形地貌制图、海洋测绘、海洋信息系统、计算机图形学与海底科学研究中具有非常重要的实际应用价值。
附图说明
图1是本发明实施例中的工作流程图。图2是本发明实施例中的模型分块图。图3是本发明实施例中的原始数据查询图。图4是本发明实施例中的图层综合制作多波束水深图。
具体实施方式
下面结合附图和实施例对本发明的做进一步说明。实施例11、一种基于DDM特征提取的多波束水深图构建方法,包括下列步骤:本发明实例完成的技术流程图见图1;步骤1、基于原始多波束测深的数字水深模型DDM的构建使用多波束后处理系统,对多波束原始数据进行如下处理:1.1)对获取的原始多波束测深数据集合Raw={raw_i}_i=1,n进行潮位改正、吃水改正和声速改正,以消除水位和声速对多波束测深MBES的影响,形成初步处理后的数据集合Proc1={proc1_i}_i=1,n,n为自然数;1.2)建立地形趋势基面,采用自动滤波方法对Proc1={proc1_i}_i=1,n进行处理,消除粗差对MBES的影响,形成后处理数据集合Proc2={proc2_i}_i=1,n;1.3)采用三视图投影方式,处理Proc2={proc2_i}_i=1,n,进一步消除多种来源海洋噪声对MBES的影响,进而形成全波束离散的水深数据Proc3={x_i,y_i,z_i}_i=1,m,m为自然数;1.4)采用距离反比加权方法IDW(Inversexa0Distancexa0Weighted)方法,对Proc3={x_i,y_i,z_i}_i=1,m进行数字建模,形成数字水深模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="700"/>K为该模型行数,L为该模型列数,K和L均为自然数,x_k,l,y_k,l和dep_k,l分别为模型数据点的坐标和水深值;步骤2、基于DDM建立坡度与二阶导数复合模型2.1)采用八邻点方法,计算每个模型网格点的坡度值slp_(i,j): slp ( i , j ) = 1 8 Σ i = k k + 1 Σ j = l l + 1 a tan ( Δz Δd ) ]]>Δz=|dep_(i,j)-dep_(k,l)|,Δz为模型中编号为(i,j)和(k,l)两点的水深差值;ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="700"/>Δd为模型中编号为(i,j)和(k,l)两点间的距离;2.2)采用八邻点方法,计算每个模型点的二阶导数值sec(_i,j): sec ( i , j ) = 1 8 Σ i = k k + 1 Σ j = l l + 1 a tan ( Δδ Δd ) ]]>Δδ=|slp_(i,j)-slp_(k,l)|2.3)形成复合数字模型 DSSM = { dss i , j = ( ddm i , j , slp i , j , sec i , j ) } i = 1 , k j = 1 , L ]]>步骤3基于DDM分块以提取特征点本发明实施例中的模型分块图见图2;3.1)模型分块:以边长为d的正方形隔行错位对模型DSSM分块;可根据制图比例尺确定d=[0.018×scale,0.03×scale],其中scale为比例尺值;d也可由用户直接指定;隔行错位是指第i行以d为间距分块模型网格,至第i+Δi行时,首先空隔0.5d倍的间距,然后以d为间距分块模型网格;隔行错位,形成菱形分块效果;通过分块,获得新的菱形网格模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="700"/>每个分块dss1_I,J由子模型网格ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="276"/>组成,K1和L1为分块模型的行数和列数,均为自然数;3.2)分块模型数值计算:3.2.1)计算 dss 1 I , J = { dss i , j } i = I , I + Δi j = J , J + Δi ]]>的平均水深 dep ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi ddm ( i , j ) , ]]>平均坡度值 slp ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi slp ( i , j ) ]]>和平均二阶导数值 sec ‾ I , J = 1 Num Σ i = I I + Δi Σ j = J J + Δi sec ( i , j ) , ]]>Num为每个分块模型中的格网点数,是自然数;3.2.2)按照IDW方法获取每个分块模型的中心点水深值Δdep_I,J;3.3)分块模型凹凸性判断3.3.1)当ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="328"/>且ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="204"/>时,该分块模型曲面是下凹的;3.3.2)当ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="334"/>且ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="204"/>时,该分块模型曲面是上凸的;3.3.3)其他情况下,该分块模型难以简单缺点凹凸性;3.4)分块模型特征水深确定3.4.1)当分块模型是凹曲面时,遍历模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="503"/>获取最大水深点dss_{i_max,j_max};如有多个最大水深点时,选取坡度值slp_i,j最大且靠近中心的水深点;3.4.2)当分块模型是凸曲面时,遍历模型ntent="drawing" img-format="tif" inline="yes" orientation="portrait" wi="506"/>获取最小水深点dss_{i_min,j_min};如有多个最小水深点时,选取坡度值slp_i,j最大且靠近中心的水深点;3.4.3)其他情况时,选取模型中心点dss_{i_cen,j_cen}作为特征点;3.5)模型特征点标识当dss_i,j为特征水深点时,设置其标记为1,否则设置其标记为0;步骤4基于特征点查询原始多波束水深4.1)遍历数据集合Proc3={x_i,y_i,z_i}_i=1,m,按照每个水深点的坐标值(x_i,y_i)查询它作用的模型网格点dss_I,J;模型网格点dss_I,J的坐标为(X_I,Y_J);4.2)当 dis = ( x i - X I ) 2 + ( y j - Y J ) 2 < 0.5 × grid _ d ]]>时,且当dss_I,_J为特征标识值为1时,使用IDW方法,基于数据集合Proc3={x_i,y_i,z_i}_i=1,m,重新计算模型点的水深值dep_new_I,J;grid_d为模型最近两点的空间间距,使用0.5×grid_d重新建模目的在于使重建模型点水深最接近真实值,不直接使用原位点探测水深是为了避免偶然探测误差对制图的影响;本发明实施例中的原始数据查询见图3;步骤5基于图层综合制作多波束水深图使用成图系统或地理信息系统,按照等深线层、水深层和修饰层等三个分层来制作来多波束水深图;5.1)在等深线层中,根据成图比例尺,按照国标要求构建首曲线、辅助线、记曲线和记曲线上的水深标注值;5.2)在修饰层中,添加图名、比例尺、图例、图框、经纬网线和海岸线;5.3)在水深层中,遍历模型DSSM,当模型点dss_I,J特征标识为1时,且当水深值dep_new_I,J不为空值时,使用该点为特征点,当dep_new_I,J为空值时,使用dep_I,J作为特征点;5.4)当最水深图中水深点过于稀疏时,减小模型分块大小值d,返回步骤3;5.5)当最水深图中水深点过于密集时,增加模型分块大小值d,返回步骤3;5.6)多图层叠加后,如出现图形要素相互叠加而不可见现象时,或者图形要素距离太密集时,可根据要素疏密程度人工删除一些辅助的图形要素。步骤5.4)所述的稀疏为所制作的水深图中任意一点与其相邻点的距离D>2d;步骤5.5)所述的密集为制作的水深图中任意一点与其相邻点的距离D<0.5d;或者上述稀疏或密集根据需要由本领域普通技术人员指定。本发明实例中基于图层综合制图的水深图局部见图4。